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Abstract 

We consider the potential dominated era of Friedmann-Lemaitre-Robertson- 
Walker flat cosmological models in the framework of general Jordan frame scalar- 
tensor theories of gravity with arbitrary coupling functions, and focus upon the 
phase space of the scalar field. To study the regime suggested by the local weak 
field tests (i.e. close to the so-called limit of general relativity) we propose a nonlin- 
^sO ! ear approximation scheme, solve for the phase trajectories, and provide a complete 

classification of possible phase portraits. We argue that the topology of trajectories 
in the nonlinear approximation is representative of those of the full system, and thus 
can tell for which scalar-tensor models general relativity functions as an attractor. 
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h ' 1 Introduction 

The unknown source of observed present day acceleration of the Universe, called dark 
energy, is inspiring thorough investigations of different extensions of general relativity 
(GR) and ACDM cosmology (for recent reviews see Refs. [I]). The scalar-tensor theory 
of gravity (STG) [2] offers one such consistent possibility. Besides the usual spacetime 
metric tensor it employs a scalar field playing the role of a variable gravitational 
"constant" , to describe the gravitational interaction. In the Jordan frame STG is specified 
by two functions [3], e.g. a coupling u(^f) and a scalar potential V(^f). In fact, a wide 
class of theories of gravitation, including higher order theories [3], theories of variable 
speed of light [5], as well as low energy approximations of brane world models and string 
theories [6] can be cast into the general form of STG. 
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The weak field tests [7] pose a restriction to all alternative models of gravity including 
STG, since the Universe around us tends to be described by the Einstein tensorial gravity 
very precisely [5] . This means that only those STG models are physically viable which in 
their late time cosmological evolution imply local consequences very close to those of GR. 
Several authors have studied how general relativity acts as an attractor for a wide class 
of STGs, such that the Solar System weak field (PPN) constraints spontaneously come 
to be satisfied at late times [HI QUI HT) . In our recent papers [TJl E3] we have proposed 
a limiting process for the scalar field which describes scalar-tensor cosmological models 
relaxing to satisfy the Solar System constraints, with an indication for which classes of 
STGs the attractor behavior is realized. 

The methods of dynamical systems have proved to be a useful tool when explicit 
analytic solutions are hard to find. In STG cosmology several authors have performed 
the analysis for different specific choices of the coupling and potential [14], while Refs. 
[T5| IT6| [T7J [12] study the phase space and dynamics in the general case. 

The STG phase space point corresponding to the limit of GR is peculiar in the sense 
that the standard linearization process there is hampered by ratios which turn out to 
be indeterminate. In our previous studies [TJl H3] we assumed that these indetermi- 
nate terms vanish, which allows to treat this singular point as a standard fixed point. 
In the present paper we focus upon the potential dominated era of Friedmann-Lemaitre- 
Robertson- Walker (FLRW) flat cosmological models in the phase space of the (decoupled) 
scalar field \l/ and its (cosmological) time derivative \j? = II, and propose an approxima- 
tion which takes into account all possible finite values of these indeterminate ratios, thus 
preserving the leading nonlinear term in the field equations. We give a comprehensive 
description of the phase space trajectories of the approximate nonlinear system near the 
GR point and classify all trajectories allowed by the parameters of the theory. While the 
topology of trajectories differs in linear and nonlinear approximation, there is a corre- 
spondence in the final asymptotics, i.e. whether the trajectories end up at the GR point 
or are repelled from it. We argue that the nonlinear system accurately captures the key 
properties of the full system of STG equations near this singular point, and the topology 
of trajectories of the nonlinear approximation is representative of those of the full system. 
Therefore for any given STG model with a reasonable coupling cj(\E') and scalar potential 
V(ty) our results predict whether GR is an attractor. 

The paper is organized as follows. In the next section we review very briefly the 
scalar-tensor theory of gravity in the potential dominated era, write down the field equa- 
tions of FLRW cosmology in the form of a dynamical system, and make some general 
remarks about the phase space including the singular GR point. In Sect. [3] we introduce 
the approximation method and present linear and nonlinear systems of equations in the 
neighbourhood of the singular point. In Section [4J we present solutions (phase trajecto- 
ries) of the nonlinear system and a classification of trajectories, summarized in Table 1 
and illustrated on Figure 1. In Section [5] our claims are backed up by a simple example. 
Finally, Section [6] provides a summary and some remarks for future work. 



2 



2 Full equations and the phase space 

We consider a general scalar-tensor theory in the Jordan frame given by the action func- 
tional 



S = — - / d^xy/^g 



2k 



VR(g) - ^p-V^Vp^ - 2« 2 V(tf) 



(We have not included the matter contribution to the action, i.e. we consider only the 
potential dominated epoch.) Here uj(^f) is a coupling function (we assume 2u(^f) + 3 > 
to avoid ghosts in the Einstein frame, see e.g. Ref. [H]), V($) > is a scalar potential, 
V M denotes the covariant derivative with respect to the metric g^, and k 2 is the non- 
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variable part of the effective gravitational constant In order to keep it positive we 
assume that < $ < oo. 

The field equations for the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) line 
element 

ds 2 = -dt 2 + a{tf (dr 2 + r 2 (d6 2 + sin 2 9d<p 2 j) (2) 

read 



# 6* 2 ^ ' 3 * 



\l> 1 \[> 2 \1> K 2 

2H + 3H 2 = -2H---—u^)-- + -V(V), (4) 



■qj 2W * ^ 

•qj = -3Hty -— -V^ * 1 



2w(^) + 3 d# 2w(^) + 3 



(5) 



where H = a/a. 

By defining = II, the equations (E])-© can be considered as an autonomous dy- 
namical system, which is characterized by three variables (q/,U,H). However, one of 
them is algebraically related to the others via the Friedmann equation ([3]). We choose to 
eliminate the Hubble parameter H and thus reduce the system to two dimensions (\l/, n). 
Upon introducing the functions 

the dynamical system reads 



* = n, (7) 

(8) 



H = ( A + ^)(2^) + 3))H 2 



±2^V3(2- (*) + 3)H 2 + 12«W(*) H + ^ + g . 
Its regular phase trajectories and fixed points have been considered in Refs. [T2J IT7] . 



3 



Even if the scalar field rules the cosmological evolution, in the context of Solar Sys- 
tem experiments we can reasonably assume that the energy density of the potential is 
negligible in comparison with the local matter density. Then the standard PPN analysis 
gives a condition for the present cosmological background value of the scalar field [TJ, 
(a) 2uj(1)+3 ~~ ^ 0- ft ^ s a ^ so evident that the effective gravitational constant is virtually 
immutable in time [10J, which provides the second condition, (b) = II — > 0. Let us 
denote ^ the value of the scalar field where the coupling function u>(^/) has a singular 
"peak", or more precisely 2aj( -J~ ) +3 = 0. Also, let II* be its vanishing time derivative, 
II* = 0. Then we may cautiously call the phase space point = \I/*, II = II*) a 'GR 
point', since the STG solutions which can pass the local weak field tests and behave close 
enough to the ones of general relativity in terms of local observations, necessarily lie in 
the vicinity of this point. 

The value \I/* poses a caveat, though. On the right hand side of Eq. the (2u;(\l/*) + 
3)II 2 terms are diverging if II ^ and indeterminate if II = II* = 0. This means that 
the whole set \l/ = \I/* at arbitrary II is excluded from the (open) domain of definition of 
Eq. ([8]). However, in what follows we find solutions of general approximated equations 
(Sect. HI) and of a specific full equation (Sect. [5]) for phase trajectories which smoothly 
reach or pass through the point (\l/*, II*). Thus we are justified to add the GR point to 
the open domain of definition as a boundary point. (Note that analogous conclusion can 
be inferred from our analysis of regions in the phase space accessible to phase trajectories 
in a general STG [12].) 

The phase portrait of the dynamical system ([7]), ([8]) is drawn by the solutions (trajec- 
tories) of 

>m / A + I^)(2 W (^) + 3)")n (9) 



±\ J3(2u(V) + 3)IP + 12k 2 W(#) + 11 ' * ' 



2^ v v K ' ' K ' (2u;(tf) + 3)11 ' 

while the direction of the flow along them is determined by Eq. ((7j as usual (towards 
increasing \I/ for II > and towards decreasing \1/ for II < 0). Realizing that gives 
the slope of the tangent to a phase trajectory, a few characteristic features of the phase 
portrait outside the singular point (\l/*, II*) can be immediately inferred. 

First, on the horizontal line II = (\l/-axis) the tangents of trajectories are vertically 
aligned if V4 7 ( 1 I / ) ^ 0, since ^ diverges due to the last term in Eq. ([9]). The inverted 
derivative ^-|n=o vanishes, while the sign of ^|| n=0 ~ indicates the direction of the 
flow along the trajectories: passing from II > to II < if ^||n=o < and vice versa if 
^g||n=o > 0- The case of quadratic potential is special, as now W / ( , I / ) = 0, and Eqs. (J7j), 
OH]) reveal that besides the singular (indeterminate) point (\1/*,II*), all points on the line 
II = are fixed points, i.e. the trajectories do not pass through, but either begin or end 
there. 

Second, approaching on the vertical line \1/ = \1/* the tangents of phase trajectories 
turn again vertical due to (2w(\l/) + 3) blowing up as \1/ — > \I / *. As has been argued 
above, the line ^ = ^ (II ^ 0) itself does not belong to the domain of the definition of 
the system, and here we see that the trajectories acknowledge this fact by not running 
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into that line, but radically deflecting "up" or "down" instead, depending on the sign of 

dU I 

These two features qualitatively control the behavior of trajectories around the GR 
point. Consider for instance the region ^ < T'*, IT > where the flow is directed towards 
First, if ^Jr| n =o < at least some of these trajectories will cross over to the II < 
belt and then flow away from \l/*, effectively displaying a "saddle" type of behavior. On 
the other hand, if ^j||n=o > there is no other option for the trajectories but to persist 
in flowing towards \I/*, while additional trajectories coming over from the II < region 

join their course. Second, near the \1/ — \I/* line if -r^- 1 <j/ >vi/^ is positive, the flow is pushed 

"upwards" to II — > oo, while if ^l^-^ is negative, the push is "downwards" towards 
smaller values of II. In the latter case the outcome again depends on how the flow is 
directed on the II = line, viz. if ^p-| n=0 < the trajectories can cross over to the II < 
region at any ^ ^ ^ and flow away, while if ^§| n=0 > the trajectories would have 
nowhere else to go but to hit the point U*). A similar reasoning can be put forth 
for \l/ > as well. It may also happen that as the potential V(^) varies, the quantity 
|jj||n=Q ~ w\v) ma y ^ e P os itive or negative depending on the value of $ and the global 
picture gets rather complicated. 

What can be said about the trajectories at the singular (indeterminate) point (\I/*, 
H*)? The logic of the phase space tells that passing through this point is possible from 
a region of the phase space where the flow is directed towards the point into a region 
of the phase space where the flow is directed away from the point. Under a reasonable 
assumption that (at least) the (physically relevant) solutions are continuous and smooth 
there are two such possibilities. First, the solutions may "slip through" ^ from the 
region # < IT > (* > II < 0) to the region * > II > (^ < II < 0) if 
= 0, i.e. the tangent of the trajectory is aligned horizontally. Second, the 
solutions may "bounce back" from \I/* from the region ^ < T^, IT > > II < 0) to 
the region * < II < > n > 0) if f| | ^-^n-s-ii* = 0, i.e. the tangent of the 
trajectory is aligned vertically. Any trajectory hitting this point under a tangent which is 
neither horizontal nor vertical can not pass through the point, but must terminate there. 



3 Approximate equations 

The equations ([7]), ([8]) cannot be integrated without specifying the two arbitrary functions 
u(^>) and V(ty). But being interested in the behavior of solutions close to the GR point 
(\I/*, U*) we can still proceed by considerning an approximation which maintains the key 
properties of the full system near this point. Although the full equations become singular 
(indeterminate) at (\l/*, II*), we assume that the Taylor expansions of the functions oo(^f), 
V(ty) are possible there. 

Let us focus around the GR point, 

V = V ir + x, U = U, + y = y, (10) 

where x and y span the neighbourhood of first order small distance from (\I/*, II*). As 
phase space variables x and y are independent from each other, and so their ratio y/x 
is indeterminate at (x — 0, y — 0). The meaning of this indeterminacy is perhaps 
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better illuminated in the polar coordinates (p, 6), where the radius p is a first order small 
quantity, but y/x = tan 9 G (—00, 00) becomes infinitely multivalued at the origin. 
We can Taylor expand 

1 1 

- + A k x + ...f*iA ic x, (11) 



2u(^) + 3 2w(**) 
and 

^ + 3)112 -ot£^. = £ (i+o(x)) -£' « i2) 

where A+ = t4(\E^). In order to keep the expansion under better control we have intro- 
duced here two additional conditions: (c)A* 7^ 0, and (d) 2 uj(i)+3 * s different iable at 
(A* and higher derivatives do not diverge) [T2l IT3] . Although these assumptions some- 
what constrain the possible forms of u, we are still dealing with a wide and relevant 
class of theories. In fact, the set (a)-(d) guarantees, that the second condition for the 
cosmological background value of the scalar field arising in the PPN analysis of STG [7], 



— > 0, is automatically satisfied [12J. 



(2w(*)+3) 3 d* 

To simplify the notation, let us denote the values of some functions at (\&*, IE) as 



Cl . ± ^™ c, S A.W„ (13) 

where W± = W(^+) and V(\&*) > 0. The three constants A+, W±, C\ determine the 
leading terms in expansions of the two functions oj(^>), V(^) which specify a STG. Now 
the expansion of the solution for H of the Friedmann constraint reads 




H = -^±^^ + 3)^ + ^0 (14) 



3k 2 d fV\ 3 y 2 

^ x + o T , — + 



This explains the introduction of the ± sign in the definition of C\ in Eq. ( 1T3|) . as near 
the GR point (x = 0, y = 0) a positive constant, Ci > 0, describes an expanding de Sitter 
Universe, while a negative one, C\ < 0, describes a contracting de Sitter Universe. 

Having outlined the method of approximation in the neighbourhood of II*), let 
us apply it for the system (J7D, dHJ). If we assume, motivated by the condition (b), that 
physically relevant trajectories linger in the region close to the x-axis, i.e a = tan# being 
first order small, then in the expansion of Eq. (JSJ) only the terms linear in x and y survive 

2 _____ 

at the first order, while terms like (cf. Eq. f )12p ) can be dropped. This was the 
assumption implicit in our earlier analysis [121 US] ■ m this case, denoting the variables x 
and y, the approximation of (JTJ) , (JHJ) yields a linear system 

x = y, (15) 
y = C 2 x -Cxy, (16) 

which, of course, is equivalent to a general second order linear homogeneous differential 
equation 

'x + d x-C 2 x = 0. (17) 
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Its phase space analysis is well known, there is a fixed point at (x = 0,y = 0) whose type 
depends on the values of the constants C\ and C 2 (see e.g [T9"]). 

2 

However, in a more general case we must recognize the term as being the same 
order as x and y. In other words, we consider all finite values of tan 9, and exclude only 
its infinite value on the y-axis which is outside the domain of definition of the system 
as said before. Thus, keeping the term ^- in the approximation of (j2J), (jHJ), we obtain a 
nonlinear system 



x 



y 



2x 



C\ y + C 2 x . 



The corresponding second order nonlinear differential equation reads 

;>2 



x + C\ x — C 2 x 



X 

2x 



(18) 
(19) 

(20) 



Note that as distinct from eqs. ( 1151) . ( |T6i) . in this case point (x = 0, y = 0) is not a fixed 
point with x = 0, y = any more, but it is a singular point with an indeterminate and 

2 

possibly multivalued term |- in eq. ( TT9|) . 



4 Phase trajectories 

The phase trajectories for the nonlinear approximate system (fT8"j) . ( JT9"j) are determined by 
the equation 

(21) 
(22) 



dy 

dx 



y_ 

2x 



x 



Cx + -C 2 . 

y 



Its solutions 



\x\K 



-y 

T 



1 1 C\yx — C 2 x 2 



exp(-Ci/(u)) , u = - 



depend on the sign of the expression C\ + 2C 2 = C, as the function f(u) is given by 



/(«) 



C 



In 



« + Ci - v 7 ^ 



u + Ci + >/C 



u + Ci 

2 r. 



c\ V 



arctan^ + nvr 
CI 



if C>0, 

if c = o, 

if C<0. 



(23) 



Here is a constant of integration which identifies the trajectory according to initial 
data (xq, yo). Note that if we choose initial conditions from the allowed region (tan# is 

2 

finite), then our premise — ~ y at deriving approximate equations ( ITS]) . (lT9"j) is always 
valid, i.e. we get a small constant of integration, K < 1. 

In general, the right hand side of Eq. (1211) can be written as a quotient of two second 
order homogeneous polynomials; a qualitative classification of the solutions of differential 
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equations of this type was given by Lyagina [20] long time ago. In a nutshell, the phase 
portraits for different values of the constants C\ and C2 classify according to the number of 
sectors which form on the phase space around the origin (x = 0, y = 0), and the topology 
of trajectories which inhabit these sectors. The sectors are separated by the boundary 
x = and invariant directions. The latter are lines y = kx where the constant A; is a real 
solution of an algebraic equation 

k = k 2~ c ' + T' (24) 

i.e straight trajectories y = (—C% ± \/~C)x satisfying ( )2~Tj) . All possible options are listed 
in Table 1 and graphically depicted on Figure 1. 

If C > and C 2 7^ three directions divide the phase space into six topologically 
distinct sectors. The sectors can be elliptic where all trajectories start from the origin 
and get back to the origin, hyperbolic where all trajectories flow towards the origin but 
turn back before reaching it, or parabolic where all trajectories either start from afar and 
flow to the origin (stable case) or start from the origin and flow away (unstable case). 

If C\ 7^ and C 2 = 0, i.e. the potential has a special form V{^>^) 7^ 0, 

= 0, then it follows from the original dynamical system ([TS]), 
f|T9|) that the entire x-axis (the point x = excluded) is populated by degenerate fixed 
points. There are four sectors. The two sectors which contain the x-axis are special, and 
do not properly belong to neither elliptic, hyperbolic or parabolic class, as the flow there 
is dominated by the cohort of fixed points lying on the x-axis. Let us provisionally call 
these 'sectors of degenerate fixed points'. 

If C = there are four sectors. If both G\ = and C2 = 0, the x-axis consists again 
of fixed points, while the generic trajectories are parabolas 2i^|x| = y 2 . If C < there 
are no real solutions to Eq. ( 12"T|) and all we get are two elliptic sectors on both sides of 
the y- axis. 

It is worth pointing out here, that the phase portraits of the nonlinear approximation 
display the same basic characteristic features we inferred about the solutions of the full 
system fl9]) before. First, on the horizontal axis (y = 0) the tangents of the trajectories 
are vertically aligned if C 2 7^ 0, and the direction of the flow across y — is determined 
by the sign of j^\ y=0 . If C2 = the horizontal axis is populated by fixed points. Second, 
next to the vertical axis (x = 0) the trajectories turn vertical and do not cross or intersect 
with the x = 0, y 7^ line, deemed to be outside of the domain of definition of the system. 

What happens at the origin (x = 0,y = 0) where the sectors meet, needs an extra 
consideration. Inspection of the phase portraits on Figure 1 shows that in all cases there 
are multiple trajectories (identified by different values of K) which all reach the point 
in question. Although our solutions for the trajectories, Eq. ( |22|) . are given in terms 
of the phase space variables only and do not include time as an explicit parameter, the 
considerations presented in the end of Sect. 2 allow to draw some qualitative conclusions. 
First, there was a logical possibility of trajectories "slipping through" the origin so that x 
changes its sign along a trajectory. It is evident from the phase portraits that this option 
is not realized in any of the cases, as none of the trajectories has tan# = at this point. 
On the other hand, the second possibility, where the trajectories could "bounce back" 
from the origin so that y changes its sign along a trajectory, is common to all cases, for 
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there is always a class of trajectories whose tangent is vertically aligned at this point. 
Despite the fact that there seems to be loss of predictability here (the initial condition 
Xq = 0, y = does not fix the constant K uniquely), it would be natural to continue 
all such trajectories through this point keeping the same K along them. Finally, those 
trajectories which reach the origin under finite tan 9 must either begin or end their flow 
at this point, like it happens at a regular fixed point. 

It is also instructive to compare the phase portraits in the linear and nonlinear ap- 
proximations. In the linear approximation (I15p . (fTB"j) the phase diagram is determined 
by 

^ = -C 1 + *C 2 . (25) 
ax y 

Its solutions are 

K=\f + C,xy - C 2 x 2 \exp(-^-f(u)) , u = ^ (26) 

2 x 

where f(u) is given by the same expression ( |23l) . while instead of C the solution is set by 
the sign of the expression C = C\ + 4CV Now the invariant sectors are separated by the 
lines y = kx, where the constant k is a real solution of 

fc = -C 1 + S, (27) 

i.e. k = §(— Ci ± vC). There is no indeterminacy in Eqs. ( ITS"]) . (1TB]) and the point 
(x = 0, y = 0) figures as a regular fixed point. The corresponding phase portraits in the 
neighbourhood of the origin can be classified according to a standard analysis [121 113] - 
the results being presented in Table 1. 

The phase portraits in the nonlinear and linear approximation look markedly different, 
as in the linear case the line x = is not a boundary, and the trajectories for nodes, focuses 
etc. are not hindered from crossing it. Yet, modulo a factor of 2 in C vs 4 in C, there is 
an overall correspondence in the classification, i.e. each distinct class in the linear case is 
matched with a distinct class in the nonlinear case (cf. Table 1). 

This correspondence in the classification occurs due to an (accidental) property that 
the nonlinear system (fTSj) . ( !T9|) for (x, y) can be formally obtained from the linear system 
( ITS"]) . (|T6|) for (x,y) by a replacement C 2 — » ^ (to get C — >■ C) and a simple transfor- 



mation x = J\x\, y = — %=. The latter can be understood as a mapping between the 

V 2y/\x\ 

corresponding phase spaces. In particular, it squeezes the whole y-axis to a single point 
(x = 0, y = 0), which manifests as the indeterminacy of the nonlinear system at this 
point. Another characteristic property of the map is the fact that the whole phase space 
(x, y) is mapped only on one half of the phase space (x,y), but covering it twice. Since 
the transformation contains absolute value \x\, there are two separate images symmetric 
with respect to the y-axis. 

Although the mapping between the linear and nonlinear systems seems to be only 
a mathematical coincidence and not a consequence of the fact that both the linear and 
nonlinear systems originate as approximations to the full system of equations (the nonlin- 
ear being a more general and refined one in this respect), we can nevertheless utilize this 
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correspondece to unravel some useful information. Namely, the trajectories which cross 
the x = line at arbitrary y in the linear case are mapped onto the trajectories which 
vertically hit the point (x = 0, y = 0), and thus, as the former pass through the x = 
line the mapping suggests that also the latter must pass through the point (x — 0, y — 0), 
hence supporting our reckoning above. Moreover, for C < or C < the trajectories 
depend not only on the real constant K or K, but also on the integer n due to the period- 
icity of arctan in Eq (1251) . Here in the linear case we have either a stable focus (spiralling 
trajectories flowing inwards into the origin) or unstable focus (spiralling trajectories flow- 
ing outwards from the origin), while n decreases by 1 on each occasion when a trajectory 
flows through the y axis. The mapping tells now, that also in the nonlinear case n must 
decrease by 1 on each occasion when a trajectory flows through the point {x = 0, y — 0), 
thus the overall picture being trajectories looping closer and closer or farther and farther 
from the origin. (On Figure 1 the diagrams 3. a and 3.c depict a trajectory with some 
fixed K but different values of n, the direction of the flow is indicated by the increasing 
number of arrows on a trajectory loop.) 

Thus (the mismatch of factor 2 vs 4 notwithstanding) besides the correspondence in 
the classification of the phase portraits, the linear and nonlinear approximations also 
share a qualitative correspondence in final asymptotic state of the flow, i.e. whether it 
ends up at the origin, or departs away from it. To summarize the results, it turns out 
that the GR point is an attractor for the asymptotic flow of all trajectories only if C\ > 
and C 2 < ( CcLSGS X.C ; 2.cl ; 3.cL ). If C\ > and C2 = all trajectories flow to the line 
^ 7^ II = instead (case l.b). If C\ = and C2 < all trajectories loop through the 
GR point oscillating back and forth (nonlinear case 3.b), or if C\ < and C2 < — 2 they 
oscillate further and further (nonlinear case 3.c). For the rest of the values of C\ and C2 
all trajectories eventually flow away from the GR point. 

5 A special example 

To further illustrate how well the full system and its approximations near the GR point 
fit together, let us consider a STG with a specific coupling function and potential, 



From the condition 2lo + 3 > it follows that \P < 1. The GR point is at ^ = 1, IT = 0, 
while the rest of the border line, ^ = 1, II* ^ 0, is deemed to be outside of the definition 
of the system. 

The full equation for phase trajectories, (jSJ), reads now 



3^ 



V(tf) = 0. 



(28) 



2(1 - *) 



dli 



3_n 

2l> 



(-1± 



2(1 -tf) 




(29) 



Fortunately it is amenable to integration, the solution is 



IT| = KWl - #(1 ± Vl - ^) 3 , 



(30) 



where K\ is a constant of integration. 
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No. Parameters Fixed point, Topology of trajectories, 

LS:N = 4; NLS:iV = 2 linear system nonlinear system 



C\ + NC 2 > 



La Ci > 
Lb d > 



C 2 > saddle 2 hyperb., 2 st. & 2 unst. parab. sectors 

C 2 = non-hyperbolic 1 stable & 1 unstable parabolic sector, 

2 stable sectors of degenerate fixed points 



l.c 


Ci 


r* 2 

^ u AT 


< C 2 < 


stable node 


2 elliptic, 4 stable parabolic sectors 


l.d 


Ci 


= 


c 2 > 


saddle 


2 hyperb., 2 st. & 2 unst. parab. sectors 


l.e 


Ci 


< 


c 2 > 


saddle 


2 hyperb., 2 st. & 2 unst. parab. sectors 


l.f 


Ci 


< 


c 2 = o 


non-hyperbolic 


1 stable & 1 unstable parabolic sector, 



Lg c x < 



4 < c 2 < 



iV 



unstable node 



2 unst. sectors of degenerate fixed points 
2 elliptic, 4 unstable parabolic sectors 



2. C^ 2 + NC 2 = 



2.a Ci > 
2.b d = 

2.c d < 



d = - 



c 2 
c 2 







stable node 2 elliptic, 2 stable parabolic sectors 
free motion 2 stable & 2 unstable parabolic sectors 



c? 



— unstable node 2 elliptic, 2 unstable parabolic sectors 



3. C\ + NC 2 < 



3.a Ci > 
3.b Ci = 

3.c Ci < 



C 2 < - 



^ 2 ^ — stable focus 2 elliptic sectors 
C 2 < centre 2 elliptic sectors 

C 2 < — -± unstable focus 2 elliptic sectors 



Table 1: Types of fixed points for linear system (LS) and the topology of trajectories for 
nonlinear system (NLS). Definitions: C = C\ + NC 2 , where iV = 4 for linear system and 
N = 2 for nonlinear system. 



11 




12 



Zooming to the vicinity of the GR point by (fTUj) . 1 — \I> = x > 0, II = y, the solu- 
tion (13"0]) approximates to a parabola \y\ = Kiy/x, which by identifying the integration 
constants K\ = 2K is in perfect accord with the solution y 2 = 2Kx of the nonlinear ap- 
proximate equation (121]) with C\ — 0, C2 = (case 2.b). This provides further evidence 
that the behavior of the full dynamical system ([7]), (JSJ) near the GR point is adequately 
approximated by the nonlinear system ffl8|) . ffT9|) . Note also, that in the linear approxima- 
tion we get a free motion, i.e. straight lines y = \/K } which qualitatively gives a correct 
asymptotic state for the trajectories (flow away from the GR point), but clearly does not 
provide a faithful phase portrait in comparison with the full solution. 



6 Summary and Discussion 

This paper considers general scalar-tensor gravity (STG) in the Jordan frame with arbi- 
trary coupling u>(^f) and potential V(^f). We have presented and justified an approximate 
theory for the behavior of the scalar field in flat FLRW cosmological STG models in a 
potential dominated era and in the regime where the local weak field experiments are 
satisfied. In terms of the phase space (\l/, \1/ = IT) the latter is understood as the neigh- 
bourhood of the 'GR point' (\1/*, II*), defined by (a) 2uJ (^ i ,)+3 = ^> O 3 ) = 0. We propose 
that if (c) ^( 2cJ (fr ) + 3 ) 7^ an d (d) the higher derivatives of 2uj ^ do not diverge, then 
in the neighbourhood of the GR point the nonlinear system (JIB]) . (|T9|) can be considered 
as an adequate approximate description of the full dynamical system ^ , (jSJ) , since both 
are endowed with the same characteristic features. The phase portraits, summarized in 
Table 1 and depicted on Figure 1, typically show many trajectories passing through the 
GR point either once on multiple times. In the end, only if 

does the GR point function as an asymptotic attractor for the flow of all trajectories in 
the vicinity. 

These analytic results could not have been predicted by numerical simulations, as the 
numerical calculations become rather problematic near the GR point due to the indeter- 
minacy present in the equations. It would be very intersting to study how the different 
looping behaviors through the GR point manifest themselves in terms of observational 
predictions. It would also be of obvious physical relevance to extend the analysis to the 
matter dominated case, although the treatment of the problem would face a difficulty of 
having an additional phase space dimension to deal with. 
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